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I.  INTRODUCTION 


The  research  described  in  this  report  is  one  of  the  first  attempts  to 
construct  a  model  for  the  temporal  propagation  of  conductivity  enhancements 
associated  with  the  auroral  ionosphere  during  substorm  breakup.  To  place  the 
research  into  a  proper  perspective  we  shall,  in  Section  III,  first  review  some 
of  the  observational  facts  relating  to  substorm  morphology.  Following  that,  we 
shall  discuss  our  results  and  their  relation  to  observations. 


II.  DESCRIPTION  OF  RESEARCH 

The  substorm  phenomenon  is  a  global  process  exhibiting  complex  dynamical 
features.  It  intrinsically  involves  the  interaction  of  coupling  between  the 
solar  wind,  magnetosphere,  and  ionosphere.  Phenomenoligically,  the  evolution  of 
most  substorms  can  be  broken  down  into  a  series  of  distinct  phases.  Akasofu 
(1977)  has  schematically  characterized  the  growth  and  decay  of  auroral  substorm 
as  shown  in  Figure  II.  1.  The  time  between  T=0  and  5  minutes  denotes  the  onset 


of  a  substorm.  It  is  characterized  by  an  initial  brightening  and  poleward 
expansion  of  a  quite  auroral  arc  (typically  the  most  equatorward  form).  Indeed, 
this  poleward  motion  of  the  auroral  arc  is  often  described  as  a  "poleward  leap" 
since  the  process  has  been  observed  to  occur  with  speeds  of  1  to  10  km/sec. 
Associated  with  the  substorm  is  the  generation  of  Pi2  pulsations.  Following  the 
initial  or  growth  phase  of  the  substorm,  the  phenomenon  evolves  in  the  manner 
shown  in  Figure  1C  and  ID.  This  complex  "breakup"  phase  of  the  substorm  is 
characterized  by  the  formation  of  a  feature  called  the  Westward  Travelling  Surge 
(WTS).  As  shown  in  Figure  II.  1,  within  30  minutes  after  onset  the  auroral 
structure  exhibits  a  twisting  and  simultaneous  propagation  towards  the  west.  We 
shall  discuss  the  properties  of  the  WTS  in  more  detail  below  when  we  describe 
our  efforts  to  model  the  phenomenon.  Following  the  breakup  phase,  the  substorm 
enters  its  final  or  "recovery  phase"  wherein  the  auroral  system  moves  back 
toward  a  configuration  characteristic  of  quiescent  conditions. 

We  can  now  use  the  brief  "definition"  of  the  substorm  given  above  to  put  our 
research  into  a  coherent  framework.  It  was  only  within  the  last  several  years 
that  systematic  morphological  properites  were  associated  with  WTS.  Several 
important  characteristic  features  of  the  surge  are: 

a)  The  structures  can  propagate  in  the  ionosphere  at  speeds  of  1  to  30 
km/sec  (Opgenoorth  et  al.,  1983:  Yahnin  et  a!.,  1983;  G.  Rostoker, 
private  communication.).  Indeed,  these  speeds  are  typically  an  order 
of  magnitude  higher  than  E  x  B  drift  velocities  and  are  moreover 
often  in  the  opposite  direction.  The  instantaneous  speeds  of  30 


km/sec  have  been  observed  in  the  leading  branch  of  the  westward 
electro  jet.  This  "horn  region"  is  westward  of  the  WTS  fold  of  Figure 
1C. 

b)  The  propagation  of  the  WTS  moves  in  steps  rather  than  continuously 
(Wiens  and  Rostoker,  1975;  Pytte  et  al.,  1976). 

c)  Correlated  with  the  WTS  is  the  appearance  of  energetic  precipitating 
electrons  (typically  1  to  10  Kev).  These  fluxes  pervade  the  high 
conductivity  region  associated  with  the  surge  (Meng  et  al.,  1978). 

d)  At  the  end  of  the  surge  there  exists  an  intense  upward  field  aligned 
current  carried  by  energetic  precipitating  electrons.  This  current 

is  concentrated  within  an  area  of  roughly  100  km  x  100  km  with  an 
average  intensity  of  1-10  yAmp/m^.  This  feature  can  be  modeled  as  a 
line  current  carrying  approximately  10^  Amps  (Inhester  et  al.,  1981). 

It  should  be  noted  that  the  head  of  the  surge  corresponds  to  the  fold 
that  develops  in  Figure  1C  and  propagates  as  shown  in  Figures  ID  and 
IE. 

e)  There  is  clearly  an  observed  relationship  between  the  WTS  and  the 
generation  of  Pi2  pulsations  during  substorm  onsets  (Rostoker  and 
Samson,  1981;  Samson  and  Rostoker,  1984). 
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The  foundation  for  research  proposed  here  is  contained  in  two  papers  which 
describe  a  model  for  the  propagation  of  the  WTS  and  the  associated  generation  of 
Pi2  pulsations.  These  works  are  included  herein  as  Appendices  A  and  B.  Both 
are  the  results  of  a  collaboration  between  myself,  P.  Rothwell  and  L.  Block. 

For  reference  we  will  denote  as  Paper  1  the  work  in  Appendix  A  entitled,  "A 
Model  for  the  Propagation  of  the  Westward  Travelling  Surge"  (JGR  89,  A10,  8941, 
Oct.  1984).  Paper  2,  entitled  "Pi2  Pulsations  and  Westward  Travelling  Surge" 
(Appendix  B),  has  been  accepted  for  publication  in  JRGR.  In  describing  our 
previous  research  we  will  concentrate  upon  summarizing  the  major  results  and 
refer  to  the  appendices  for  the  mathematical  details. 


Steady  state  models  of  the  WTS  have  been  developed  using  preliminary  ground- 
based  observations  (Hughes  and  Rostoker,  1979;  Rostoker  and  Hughes,  1979;  Tighe 
and  Rostoker,  1981).  A  comprehensive  static  model  of  the  electrodynamic 
structure  within  the  WTS  was  developed  by  Inhester  et  al.,  (1981).  This  model 
is  described  in  detail  in  Section  2  of  Paper  1.  Essentially,  one  can  model  the 
WTS  at  a  given  instant  of  time  as  a  Cowling  conductivity  channel  in  the 
ionosphere.  The  channel  is  described  as  a  slab  aligned  in  the  "east-west" 
direction.  Examination  of  Figure  II. 1  above  indicates  that  the  slab 
approximation  is  a  reasonable  one  for  certain  aspects  of  the  complex  WTS 
phenomenon  (see  eg..  Figure  IB  and  horn  region  of  Figures  1C  and  ID).  Modeling 
the  surge  head  of  Figures  1C,  D  and  E  by  a  slab  is  clearly  only  a  first 
approximation.  An  improvement  of  this  model  is  one  of  the  research  thrusts 
proposed  here. 

Kan  et  al.,  (1984)  have  pointed  out  the  importance  of  the  degree  of 
ionospheric  current  closure  on  the  poleward  boundry  of  the  enhanced  conductivity 
region  associated  with  the  WTS.  Closure  is  governed  by  the  parameter  a,  which 
is  a  measure  of  th  degree  to  which  the  net  ionospheric  current  is  continued  into 
the  magnetosphere  via  field  aligned  currents  at  the  slab  boundries  (see  Section 
2  of  Paper  1  for  more  details).  In  our  context  strong  current  closure  (a  =  1) 
implies  full  continuation  of  the  ionospheric  current  into  the  magnetosphere. 
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In  Paper  1  we  develop  a  dynamic  model  for  the  propagation  of  the  WTS.  We 
explicitly  utilize  the  closure  concept  and  incorporate  it  with  elements  of  the 
Inhester-Baumjohann  slab  model  (Baumjohann,  1983).  The  motion  of  the  WTS  in  the 
midnight  sector  is  considered  to  be  controlled  by  three  mechanisms.  These  are: 

1)  the  energy  and  intensity  of  the  precipitating  electrons; 

2)  the  electron-ion  recombination  rate; 

3)  the  degree  of  current  closure  on  the  boundaries  of  the  slab. 

Both  the  closure  parameter  a  and  electron  precipitation  energy  represent  the 
magnetospheric  input  into  the  surge  dynamics.  At  present  these  inputs  are  not 
self-consistently  incorporated  into  our  models.  A  consideration  of  this  issue 
will  be  one  of  the  important  tasks  in  our  proposed  research  program. 

The  main  idea  of  Paper  1  can  be  expressed  as  follows  (L.  Block,  private 
communiation,  1985).  If  a  Birkeland  current  moves  (e.g.,  due  to  some 
magnetospheric  action)  the  connecting  ionospheric  current  may  encounter  some 
relatively  low  conductivity  regions.  This  could  demand  a  higher  voltage  from 
the  magnetospheric  generator,  unless  the  high  conductivity  region  is  able  to 
follow  the  motion  of  the  Birkeland  current.  The  field  aligned  Birkeland  current 
is  itself  an  ionizing  mechanism  through  the  associated  energetic  electron 
precipitation  especially  in  the  upward  current  region.  Hence,  as  the  upward 
current  moves  it  can  give  rise  to  a  synchronous  motion  of  the  conductivity 
channel  provided  it  doesn’t  move  to  rapidly.  One  can  derive  some  relations 
between  the  Birkeland  current  intensity  and  the  speed  at  which  its  footprint 


moves  in  the  ionosphere.  The  validity  of  this  approach  depends  upon  the 

assumption  that  one  can  obtain  some  relation  between  the  field  aligned  current 

density  and  the  ionizing  efficiency  of  the  associated  precipitating  electrons. 

This  was  in  fact  done  in  Paper  1  and  a  wave  equation  was  derived.  The  velocity 

of  propagation  is  proportional  to  VdQH,  where  V  is  the  magnitude  of  the  total  E 

u  d 

x  B  drift  velocity  and  QH  is  the  height  integrated  ionization  efficiecy  for 
precipitating  electrons  at  the  conductivity  gradient.  It  is  a  characteristic  of 
QH  that  the  velocity  is  greatly  enhanced  when  the  average  precipitating  electron 
energy  increases  from  I  Kev  to  10  Kev.  The  variations  of  QH  with  energy  are 
shown  in  Table  1  and  Figure  3  of  Paper  1.  As  an  example  we  note  that  a  10  mV/m 
ionospheric  electric  field  implies  a  propagation  velocity  of  roughly  3.7  km/sec 
for  1  KeV  incident  electrons  with  10  Kev  electrons  will  produce  a  velocity  of  34 
km/sec.  Thus,  in  applying  this  propagation  mechanism  to  our  slab  model  of  the 
WTS  we  predict  that  jumps  in  the  surge  velocity  will  occur  whenever  there  is  a 
hardening  of  the  energy  spectrum  of  the  precipitating  electrons  associated  with 
the  upward  field  aligned  current  region  of  the  WTS. 

The  simple  dynamical  slab  model  of  Paper  1  predicts  a  wide  range  of  complex 
phenomena  associated  with  the  poleward  leap,  WTS  propagation  and  recovery  phase 
of  a  substorm.  In  this  model  the  important  parameters  are  the  ratio  of  Hall  to 
Pedersen  conductivities,  the  degree  of  ionospheric  current  closure  into  the 
magnetosphere,  and  the  energy  spectra  of  the  precipitating  electrons.  We  will 
conclude  this  subsection  by  summarizing  the  major  results  of  Paper  1.  These 


The  direction  of  the  WTS  depends  strongly  on  the  degree  of  current 
closure  on  the  poleward  boundary.  The  sensitivity  of  the  surge 
direction  to  closure  depends  on  the  ratio  (R)  of  the  Hall  to  the 
Pedersen  conductivities. 

The  magnitude  of  the  surve  velocity  is  sensitive  to  the  energy  spectrum 
of  the  precipitating  electrons  and  weakly  dependent  on  R.  The  ratio  of 
the  surge  velocity  to  the  measured  drift  velocity  is  independent  of  the 
degree  of  closure  at  both  the  surge  head  and  the  northern  boundary. 

The  initation  of  the  expanse  phase  of  the  substorm  can  be  explained  by 
assuming  that  the  initial  arc  brightening  arises  from  a  sudden  hardening 
of  the  precipitating  electron  energy  spectrum  as  its  poleward  boundary. 

Inclusion  of  electron-ion  recombination  effects  highlights  the  role  of 
the  precipitating  current  intensity  in  modulating  the  surge  propagation 
and  explains  the  equatorward  retreat  of  the  surge  during  the  substorm 
recovery  phase. 

Details  of  the  surge  propagation  depend  on  how  the  magnetosphere  and 
ionosphere  are  coupled  as  reflected  by  the  quantities  a  and  Q.  Hence 
the  energy  source  is  clearly  located  in  the  magnetosphere  and  a  complete 
description  of  substorm  phenomena  must  take  this  into  account. 


The  second  major  component  of  our  previous  research  concerns  an  investigation 
into  the  generation  of  Pi2  pulsations  and  their  correlation  with  substorm 
dynamics.  As  indicated  above,  the  results  of  this  work  are  detailed  in  Paper  2 
which  is  reproduced  here  as  Appendix  B.  As  in  the  previous  subsection  we  shall 
simply  summarize  the  important  results  of  this  research  and  refer  the  reader  to 
Paper  2  for  specific  derivations.  To  put  our  results  into  a  coherent  framework 
we  first  list  three  important  characteristics  in  the  Pi2  waveforms.  These  are: 

a)  The  appearance  of  Pi2  wavetrains  is  correlated  with  the  onset  of  a 
substorm. 

b)  The  Pi2  waveform  has  a  damped  quasi-sinusoidal  shape.  The  periods  of 
the  pulsations  fall  between  40  and  140  seconds.  Their  duration  is 
typically  on  the  order  of  several  minutes.  Figure  4  of  Paper  2  shows 
some  representative  data  for  the  Pi2  pulsations. 

c)  The  polarization  sense  of  the  waveform  changes  if  one  compares 
observations  inside  and  outside  the  regions  of  enhanced  conductivity 
characteristic  of  substorm  activity. 

Paper  2  refines  the  dynamics  of  the  simple  slab  model  developed  previously. 

To  incorporate  the  generation  of  Pi2  waves  we  assume  that  the  ionosphere- 
magnetosphere  coupling  can  be  modeled  by  a  transmission  line.  This  approach  was 


first  developed  by  Sato  (1978)  and  takes  into  account  the  propagation  of 
transverse  Alfven  waves  generated  by  perturbations  that  originate  either  in  the 
ionosphere  or  magnetosphere.  In  our  calculation  we  consider  the  effects  of  an 
ionospheric  perturbation  which  occurs  in  the  properties  of  the  Cowling  channel 
used  to  model  the  WTS  dynamic.  The  effect  of  the  transmission  line  coupling  is 
to  include  the  possibility  of  a  feedback  instability  which  can  cause  ionospheric 
disturbances  to  be  amplified.  This  instability  mechanism  was  first  introduced 
by  Atkinson  (1970)  and  later  refined  by  Sato  and  his  co-workers.  All  previous 
work,  however,  had  been  applied  to  the  properties  of  quiet  auroralforms.  To  the 
best  of  our  knowledge.  Paper  2  represents  the  first  attempt  to  utilize  the 
feedback  instability  mechanism  in  a  theory  of  substorm  dynamics.  The  major 
results  of  Paper  2  are: 

1)  The  use  of  the  feedback  instability  coupled  with  the  dynamical  results 
of  Paper  1  yields  a  theory  which  ties  the  generation  of  Pi2  pulsations 
of  substorm  onsets.  In  particular  we  have  developed  a  perturbation 
scheme  based  on  a  closed  set  of  nonlinear  partial  differential 
equations.  To  zeroth  order  our  thecry  is  exactly  that  of  the  bulk 
propagation  of  the  WTS  described  in  Paper  I.  To  first  order  we  obtain  a 
dispersion  relation  which  reflects  the  feeback  instability  mechanism. 

The  frequency  and  growth  rate  of  the  modes  are  shown  in  Figure  3  of 
Paper  2.  There  it  is  seen  that  the  frequency  range  is  characteristic  of 
the  Pi2  pulsation. 
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The  frequency  and  growth  (or  damping)  of  the  Pi2  waves  is  a  function  of 
WTS  velocity.  From  Paper  1  this  implies  that  the  properties  of  Pi2s 
depend  on  the  characteristics  of  the  precipitating  electrons  associated 
with  an  enhanced  conductivity  region. 

In  Paper  1  it  is  found  that  the  slab  propagates  only  if  the  field 
aligned  current  exceeds  a  critical  threshold  wherein  ionization 
dominates  over  recombination  effects.  For  parameters  characterizing 
substorm  conditions  the  critical  current  value  is  on  the  order  or 
5yAmp/m  When  the  current  threshold  is  superimposed  on  the  wave 
characteristics  of  Paper  2  it  is  found  that  the  Pi2  modes  are  damped  in 
agreement  with  physical  constraints.  This  result  will  be  discussed  in 


I 


more  detail  in  Section  IV.2. 
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A  unified  model  is  developed  for  the  propagation  of  the  westward  traveling  surge  (WTS)  that  can 
explain  the  diversity  in  the  observed  surge  characteristics.  We  start  with  the  Inhester-Baumjohann  model 
for  the  surge  region,  which  implicitly  includes  both  the  Hall  and  Pedersen  currents.  It  is  Found  that 
precipitating  electrons  at  the  conductivity  gradient  modify  the  gradient,  causing  it  to  propagate  as  a  wave 
front.  The  velocity  of  propagation  is  directly  dependent  on  the  ionization  efficiency  of  the  precipitating 
electrons  and  therefore  increases  dramatically  when  they  become  more  energetic  during  substorm  onsets. 
For  example,  we  predict  that  when  the  incident  electron  energy  changes  from  1  keV  to  10  keV  the  surge 
velocity  should  increase  from  2  km/s  to  34  km/s.  The  direction  of  the  surge  motion  depends  on  the 
presence  of  polarization  charges  on  the  poleward  surge  boundary.  This  is  related  to  the  efficiency  with 
which  the  poleward  ionospheric  currents  are  closed  off  into  the  magnetosphere  by  the  field-aligned 
currents.  Inclusion  of  the  electron-ion  recombination  rate  modifies  the  surge  propagation  velocity  and 
leads  to  explicit  expressions  for  the  conductivity  profile.  Sufficient  precipitation  current  is  required  to 
overcome  electron  ion  recombination  in  order  for  the  surge  to  expand.  When  the  precipitating  current  is 
less  than  this  threshold  the  WTS  retreats.  Therefore,  the  model  describes  the  ionospheric  response  to 
both  the  expansion  and  recovery  phases  of  the  magnetic  substorm 


I.  Introduction  Hughes  and  Rostoker  [1979],  Rostoker  and  Hughes  [1979], 

The  westward  traveling  surge  (WTS)  is  a  primary  signature  and  T‘9he  and  Rostoker  [1981]  using  primarily  ground-based 

of  substorm  onsets  [Rostoker  el  al.,  1980],  In  simplest  terms  observations.  By  incorporating  STARE  data,  Inhester  et  al. 

the  WTS  represents  the  westward  electrojet  as  it  expands  [1981]  have  developed  a  comprehensive  static  model  of  the 

westward  starting  near  local  midnight  during  the  subslorm  electrodynamic  structure  within  the  surge.  Kan  et  al.  [1984] 

expansion  phase  [ Tiyhe  and  Rostoker,  1981],  Figure  I  shows  a  have  pointed  out  the  importance  of  ionospheric  current  clo- 

typical  WTS  profile  as  measured  by  the  DMSP  satellite  F-6.  sure  on  the  poleward  surge  boundary.  Using  static  conduc- 

In  this  example  we  see  evidence  for  multiple  surge  structure  as  tivity  models  they  show  that  closure  is  responsible  for  the 

described  by  Rostoker  et  al.  [1980],  Note  that  the  shape  of  the  westward  intrusion  of  highly  conductive  regions.  In  the  pres- 

WTS  is  similar  to  that  of  oceans  waves  where  the  wave  crest  en*  dynamical  model  we  explicitly  show  how  closure  deter- 

corresponds  to  a  surge  head.  A  more  detailed  examination  of  mines  the  surge  direction  and,  in  contrast  to  their  conclusions, 

the  surge  reveals  several  characteristic  features  of  the  phenom-  hod  that  the  surge  speed  does  not  necessarily  depend  on  the 

ena.  These  include  the  following:  (1)  The  surge  moves  in  dis-  t‘me  rate  change  of  closure. 

crete  jumps  or  steps  [Wien  and  Rostoker,  1975;  Pylte  et  al.,  ln  this  PaPcr  we  develop  a  dynamic  model  for  the  propaga- 
1976],  (2)  Instantaneous  velocities  of  up  to  30  km/s  have  been  **on  ’he  WTS.  As  shown  in  the  next  section  we  explicitly 
observed  in  the  leading  branch  of  the  westward  electrojet  [Op-  utilize  the  closure  concept  and  incorporate  it  with  elements  of 

genoorth  et  al..  1983;  Yahnin  et  al.,  1983,  G.  Rostoker,  private  ,hc  Inhester  Baumjohann  model  [ Baumjohann ,  1983],  The 

communication,  1983).  However,  a  typical  surge  velocity  is  on  motion  of  the  WTS  in  the  midnight  sector  is  considered  to  be 

the  order  of  1.0  2.0  km/s  [Pytte  et  al.,  1976].  (3)  The  energy  of  controlled  by  three  mechanisms.  These  are  (I)  the  energy  and 

the  precipitating  electron  flux  associated  with  the  surge  region  intensity  or  the  precipitating  electrons,  (2)  the  electron  ion 

is  usually  in  the  keV  range  [Meng  et  at.,  1978],  (4)  At  the  head  recombination  rate,  and  (3)  the  degree  of  current  closure  on 

of  the  surge  there  exists  an  intense  upward  field-aligned  cur-  the  poleward  boundary  of  the  surge.  Closure  is  governed  by 

rent  carried  by  energetic  precipitating  electrons  This  current  the  parameter  a,  which  is  a  measure  of  the  degree  to  which  the 

is  concentrated  within  an  area  of  at  least  100  km  x  100  km  net  ionospheric  current  is  continued  into  the  magnetosphere 

(I  x  I  )  with  an  average  intensity  of  1-10  pA/m2  and  can  be  at  'he  poleward  boundary  (see  discussion  near  equation  (2)  for 

modeled  as  a  line  current  carrying  £  10s  A  [fnhesrer  et  at.,  more  details).  In  our  context,  strong  current  closure  (a  =  1) 

1981].  The  slab  model  developed  in  section  2  can  be  applied  implies  full  continuation  of  the  ionospheric  Hall  current 

to  both  the  surge  head  and  to  the  western  precursor  region.  system  into  the  magnetosphere  via  field-aligned  currents.  Both 

Steady-state  models  of  the  WTS  have  been  developed  by  3  a"d  the  precipitation  energy  represent  the  magnetospheric 

input  to  the  surge  dynamics.  The  mechanisms  describing  these 
This  paper  is  mu  subject  to  US.  copyright  Published  in  19X4  by  parameters  are  outside  the  scope  of  the  present  paper, 
the  American  Geophysical  Union  Upward  field-aligned  currents  in  the  WTS  are  most  intense 

Paper  number  4A0864  where  the  conductivity  gradient  is  largest.  These  currents  are 
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carried  by  precipitating  electrons  that  modify  the  conductivity 
gradient  through  enhanced  ionization.  A  wave  equation  is 
derived  for  the  propagation  of  the  conductivity  gradient.  The 
phase  velocity  is  proportional  to  Fdl QH  where  Vit  is  the  mag¬ 
nitude  of  the  total  £  x  B  drift  velocity  and  QH  is  the  height- 
integrated  ionization  efficiency  for  precipitating  electrons  at 
the  conductivity  gradient.  This  velocity  is  greatly  enhanced 
when  the  average  electron  energy  increases  from  1  keV  to  10 
keV  For  example,  a  10  mV/m  electric  field  implies  a  surge 
velocity  of  approximately  3.7  km  s  for  1-keV  incident  elec¬ 
trons,  whereas  10-keV  incident  electrons  will  drive  the  surge  at 
approximately  34  km  s  (see  Table  I  for  details).  We  see  no 
direct  connection  between  this  propagation  mechanism  and 
that  based  upon  an  ion-acoustic  wave  as  proposed  by  Kan  et 
al.  [1984] 

The  effect  is  to  produce  jumps  in  the  surge  velocity  when¬ 
ever  the  precipitating  energy  spectrum  near  the  surge  head  is 
sufficiently  hard  As  shown  below,  it  is  even  possible  to  have 
abrupt  northeastward  surge  motion  if  there  is  overclosure  at 
the  poleward  boundary  (i.e.,  a  >  1)  Thus  even  this  simplified 
model  provides  coherence  to  the  diversity  of  observed  surge 
propagation  characteristics. 

In  section  2  the  simplified  model  is  developed  for  the  WTS 
where  electron  ion  recombination  effects  are  ignored.  The 
magnitude  of  (he  surge  velocity  and  its  direction  are  derived. 
In  section  3  we  investigate  the  role  of  electron-ion  recombi¬ 
nation  effects.  We  find  that  the  recombination  rate  modifies 
the  surge  velocity,  determines  the  conductivity  profile  at  the 
boundary,  and  is  responsible  for  the  retreat  of  the  enhanced 
conductivity  region  and  the  associated  electrojet  during  the 
recovery  phase. 

2  A  Model  eor  the  Westward  Traveling  Surge 

We  start  with  the  slab  model  for  the  WTS  near  local  mid¬ 
night  as  given  by  Baumjohann  [1983]  and  as  shown  in  Figure 
2  In  this  section  we  only  consider  the  motion  of  this  slab.  The 
results  will  be  applied  to  the  WTS  in  the  discussion  section. 
For  clarity  the  stated  propagation  directions  are  defined  for 
the  WTS  occurring  in  the  northern  hemisphere.  An  effective 
westward  electric  he  1 1.  in  the  surge  region  drives  a  north¬ 
ward  Hall  current  This  electric  field  is  composed  of  the  large- 
scale  convection  field  and  that  produced  by  negative  charge 


current  .s 

,/H  ~  I„£0  (1) 

However,  the  net  ionospheric  current  reaching  the  northern 
h-vutuiarv  may  he  unequal  to  ,/H  depending  on  the  presence  of 
n  'undary  polarization  charges  as  indicated  in  Figure  2.  This 
boundary  polarization  charge  produces  a  souihward  polariza¬ 
tion  elect i  ic  field.  l.f.  as  shown  in  the  figure  This  electric  field 
ziv :*'  rise  to  an  associated  Pedersen  current.  J  -  I  t  Thus, 
the  net  ionosphere  current  reaching  the  northern  boundary  is 
■  =  x/H.  where  the  parameter  «  is 

.  measure  of  the  polarization  charge  at  the  slab  boundaries.  It 
s  also  ,i  measure  of  the  degree  to  which  the  net  ionospheric 
'en>  closed  into  the  magnetosphere  by  field-aligned  cur- 
•e-.is  lu.surrnng  no  ionospheric  current  closure  outside  the 
T:  -  an  undetermined  parameter  in  our  theory  that  re- 
i  lor.i.-;  j,.  re  magnetosphere  coupling.  The  value 
'  i  .  -losi.re  the  value  x  =  I  full  closure  and 


a  >  1  overclosure  (i.e.,  a  negative  polarization  charge  on  the 
northern  boundary  produced  by  excess  Birkeland  currents). 

We  define  a  coordinate  system  as  shown  in  Figure  2  such 
that  the  z  axis  points  perpendicular  to  both  the  current 
channel  and  the  earth's  magnetic  field.  The  x  axis  points 
parallel  to  the  current  channel.  When  applied  to  the  westward 
electrojet  the  x  coordinate  is  approximately  west  and  the  z 
coordinate  north.  For  simplicity,  in  the  following  model  we 
define  x  as  pointing  due  west  and  z  as  pointing  due  north. 

The  precipitating  Birkeland  current  is  the  divergence  of  J„ 
as  given  by 

7  it  =  -^£o*Ih)A- 

=  —  £0adXH/(7z  (2  a) 

where  x  is  the  coupling  parameter  defined  above.  A  meaning¬ 
ful  x  can  only  be  defined  where  there  is  a  conductivity  gradi¬ 
ent.  To  be  more  precise,  consider  current  conservation  across 
the  boundary. 

j„  dz  =  £  -  E0(Xx(z)l^z))  dz  (2b) 

If  a  =  const,  within  the  boundary  region  (z,  <,  z  <,  z2 ),  then 

J"  7 ii  dz  =  *-fH(-i)  (2c) 

This  is  another  way  of  expressing  the  overall  current  conti¬ 
nuity  at  the  northern  boundary.  An  x  <  1  implies  the  presence 
of  positive  polarization  charge.  We  associate  jM  with  the  flux  of 
precipitating  electrons  through  the  relation  j{t/e.  This  assumes 
that  the  precipitating  protons  are  not  important  in  the  surge 
region  [ Akasofu  et  al..  1969]  and  that  the  dominant  current 
carriers  are  energetic  (keV)  electrons  [ Meng  et  al.,  1978;  In¬ 
heater  el  al..  1981],  This  is  consistent  with  the  observed  en¬ 
hancement  of  keV  electrons  observed  during  the  passage  of 
the  WTS  [ Opgenoorth  et  al..  1983],  For  simplicity  we  shall 
ignore  possible  spatial  variations  in  both  x  and  £0.  The  ratio 
of  the  Hall  to  Pedersen  conductivity  (7?)  is  considered  to  be  a 
function  of  z  only. 

In  the  following  derivation  we  recognize  that  precipitating 
electrons  through  ionization  modify  the  conductivity.  The 
electron  flux.  jn/e,  may  consist  of  electrons  that  are  accelerated 
at  higher  altitudes  through  field-aligned  potential  differences. 
In  that  case  jw'e  is  the  flux  of  accelerated  electrons  incident  on 
the  ionosphere.  The  local  time  rate  of  change  of  the  Hall 
conductivity  is 

(Th  i'f  i  (cH  B)en  ft  (3) 

where  n  is  the  ion  density  that  is  related  to  the  precipitating 
electrons  through  the  ionization  efficiency  Q  [Rees,  1963;  Ja- 
sperse  and  Baau.  1982],  The  local  time  rate  of  change  of  the 
ion  number  density  is  given  by 

en/rl  =  Qi  e  -  n.n1  (4) 

where  here  Q  is  the  average  number  of  ions  produced  per 
incident  eleclron-m  and  H  is  the  appropriate  height  interval 
for  Q  and  n,  is  the  electron  ion  recombination  coefficient  in 
the  ionosphere.  (The  unil  ions, eleclron-m  when  multiplied  by 
the  unit  for  flux  (electron  m2-s)  becomes  ions  m'  s.  the  ion 
production  rate  per  unit  volume)  Equation  (4)  is  general  in 
that  an  average  Q  can  be  defined  for  any  incident  spectrum 
The  ionization  rate  is  then  this  average  Q  times  the  total 
incident  electron  current  providing  the  net  back  scattered  cur¬ 
rent  is  small.  For  example,  the  presence  of  a  parallel  electric 
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Fig.  I  The  aurora  as  photographed  by  the  DMSP-Fft  .satellite  on  January  IT  19X3  A  small  westward  traveling  surge 
(WTSl  is  observed  in  the  tenter  of  the  photograph  followed  by  a  significantly  larger  one  to  the  east  The  overall  shape  of  a 
surge  is  similar  to  that  of  an  ocean  wave  where  the  wave  crest  corresponds  to  the  surge  head  In  our  idealized  model  the 
surge  is  represented  as  a  slab  as  shown  in  Figure  ? 


field  will  effectively  eliminate  any  upward  electron  current 
contribution  [ Evans,  1974]  If  we  assume  that  only  the  ener¬ 
getic  component  (keV  electronsl  carry  (he  parallel  currenl 
then  the  Q  in  equation  (4)  is  identical  lo  the  Q  I  defined  by 

Rees  r  1963] 

We  initially  treat  the  limiting  case  where  /  (J  e  "  nrn!  in 
order  to  illustrate  the  physical  principles  involved  I  he  full 


equation  is  solved  below  in  section  3.  There  we  show  that  the 
simple  approximation  is  valid  if  j{[  »  3  m2  when  >  I 
keV  The  sign  in  front  of  /,,  in  (4)  indicates  that  positive  cur¬ 
rent  is  (lowing  away  from  the  earth  Combining  (2H4)  we 

obtain 


«’IH  el  =  -(QH L0i)iB?X.„ 
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Fig.  2.  Idealized  model  of  the  WTS  as  proposed  by  Inhester  et  at.  [1981],  The  total  external  electric  field  £0  drives  a 
westward  Pedersen  and  a  northward  Hall  current.  If  the  northward  Hall  current  is  not  fully  continued  by  field-aligned 
currents  into  the  magnetosphere  polarization  charges  build  up  on  the  northern  surge  boundary  producing  a  southward 
directed  electric  field  £p.  This  southward  electric  field  produces  a  Hall  current  in  the  same  direction  as  the  Pedersen 
current  from  the  original  electric  field  £0.  This  affects  the  direction  of  the  surge  motion  as  described  in  the  text.  Negative 
polarization  charge  may  also  build  up  at  the  surge  head  due  to  intense  electron  precipitation  [Inhester  et  al..  1981].  This  is 
treated  in  our  idealized  model  by  allowing  a  renormalization  of  £0.  In  our  coordinate  system,  z  points  perpendicular  to  £0 
and  x  is  parallel  to  £0,  which  approximately  corresponds  to  north  and  west. 


POLEWARD  BOUNDARY 


which  is  a  wave  equation  with  a  phase  velocity  given  by 

K  =  (6) 

and  where  I,  =  E0  B  Note  that  (5)  is  also  invariant  under  the 
Galilean  transformation  IH  =  IH(z),  where  z'  =  z  —  Vmt. 

Figure  3  shows  the  estimated  variation  of  QH  with  incident 
electron  energy  as  taken  from  Figure  2  of  Rees  [1963],  From 
this  figure  it  is  seen  that  a  1-kcV  monoenergetic  precipitating 
flux  implies  QH  z=  10  and  for  a  10  keV  flux  QH  ^  90.  Thus, 
when  the  energy  spectrum  of  the  precipitating  electrons  hard¬ 
ens  the  boundary  velocity  can  increase  by  factors  of  10.  Meng 
el  al  [1978]  have  measured  a  very  hard  precipitating  energy 
spectrum  in  the  surge  region.  We  argue  that  the  surge  velocity 
is  related  to  the  production  of  energetic  electrons  connected 
with  substorm  onsets.  According  to  this  mechanism  the  surge 
velocity  should  increase  dramatically  during  substorm  onsets 
as  has  been  observed  by  Samson  and  Rosloker  [1983]. 

We  now  consider  the  western  boundary  (i.e.,  the  surge 
head  i  The  dynamics  of  this  boundary  is  connected  to  the 
propagation  of  the  poleward  boundary  through  closure.  That 
i-  a  small  ,•  implies  a  large  £p  which  drives  a  westward  Hall 
cat  reiit  that  adds  to  t.i.  dedersen  current  driven  by  £0.  The 
total  westward  current  is  given  by 

3.  -  Ip£0  4  I „£„  (7) 


T\FHF  I  Height  Integrated  Ion  Production  IQH)  From 
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l  e  P'  *ent  ■  is1'  'ns  are  considered  to  he  isotropic,  mon- 
.  c.  .  ''carl,,  'lu  results  are  plotted  in  Figure  3  The  parameter 
px  i  i. re  I  ;  elKierith  the  poleward  Hall  current  closes  into 
'spher  K  "  : he  ratio  of  the  Hall  and  Pedersen  conduc- 
o.  t  I  -  d,.  treasured  isee  equation  1 1  5»l  electric  field  in  the 
■  .  pi'  i,  n  to' i1  surge  velocity  as  defined  in  (14). 


From  the  above  definition  of  3.  we  have 
3p  —  ( I  —  *)3h 

(8) 

Ip£p  =  (l  -x)£HE0 

Substituting  this  expression  for  £p  into  (7)  leads  to  an  x- 
dependent  Cowling  current 

3„  =  { 1  +  R2{  1  -  x)}E0Zp  (9) 

where 

R  = 

The  precipitating  current  at  the  head  of  the  surge  is 

3  ti  =  Jl  +  not 

3,;  =  [{<  +  R2(l  -a)}£0/R]«T„/i’x 

where  the  x  axis  is  in  the  westward  direction  and  R  is  con¬ 
sidered  independent  of  x.  Note  that  at  the  northern  boundary 
it  was  not  necessary  to  assume  R  constant  in  deriving  (5). 
Therefore.  £p  (equation  (81)  may  also  be  a  function  of  z,  and 
our  model  implicitly  allows  polarization  charge  along  the 
northern  boundary.  A  wave  equation  is  now  derived  for  the 
western  boundary  associated  with  the  surge  head  just  as  for 
the  northern  boundary.  The  resulting  phase  velocity  is 

K,  =  QH V4{1  +  3? 2( I  -x)|  R  (II) 

We  now  assume  that  the  surge  head  region  will  propagate 
in  a  direction  determined  by  the  vector  sum  of  and  F„  The 
direction  of  the  resultant  surge  velocity  based  on  this  idealized 
model  is  given  by 

tan  y,  =  l,  T'„  =  xR  I  +  RJ(  I  —  x)|  (12) 

where  y,  =  0  corresponds  to  due  west  motion.  We  see  that  the 
direction  of  the  surge  is  highly  dependent  on  the  degree  of 
closure  on  the  northern  surge  boundary  and  the  value  of  R  on 
the  western  surge  boundary.  For  zero  closure  (x  =  0)  the  surge 
moves  due  west  Foi  complete  closure  (x  =  I)  the  direction  is 
almost  due  north  [tan  y,  =  R  x  3;  y,  =;  72  ] 

The  sensitivity  of  y,  to  x  is  modulated  by  the  magnitude  of 
R  A  detailed  plot  of  y,  versus  x  for  various  values  of  R  is 
shown  in  Figure  4  Note  that  the  surge  direction  can  range 
from  due  west  to  northeast,  where  the  latter  condition  will 
arise  when  there  is  significant  overclosure  (x  >  I)  during 
periods  of  intense  electron  precipitation 


-20- 


Rothwell  et  al.  :  Propagation  of  the  Westward  Traveling  Surge 


OH  VS  E 


E(KaV) - - 

Fig  3.  Height-integrated  ion  production  efficiency  in  the  ionosphere  for  incident  precipitating  electrons.  More  ener¬ 
getic  electrons  produce  a  higher  ionization  density  near  the  end  of  range.  This  curve  is  for  an  isotropic,  monoenergetic 
incident  beam  and  was  taken  from  Figure  2  of  Rees  (t 96$].  Errors  up  to  a  factor  of  2  could  arise  from  erroneous 
estimation  of  the  appropriate  height  interval  from  Rees’  Figure  2.  A  reasonable  fit  to  this  curve  is  QH  %  5  E‘  1  ions/(inci- 
dent  electron). 


direction  of  the  westward  traveling  surge  vs 

CLOSURE  ON  THE  NORTHERN  BOUNDARY 


Fig  4  The  direction  of  the  surge  motion  as  a  function  of  the  closure  parameter  a,  which  represents  the  degree  of  Hall 
current  closure  on  the  northern  surge  boundary  The  value  a  =  I  implies  Ep  -  0  Also  a  larger  value  of  a  implies  a  more 
intense  precipitating  current  with  a  >  I  corresponding  to  a  negative  rather  than  positive  charge  buildup  on  the  northern 
boundary  The  sensitivity  of  the  surge  direction  to  a  is  dependent  on  the  ratio  of  the  Hall  to  Pedersen  conductivities  at  the 
western  surge  boundary,  which  we  denote  by  R  Note  that  yt  ~  45  corresponds  to  a  a  0  8-0.9. 
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The  magnitude  of  the  total  slab  velocity,  K„,  is  found  from 
V,,2  =  (,2  +  K.*.  Clearly,  the  resultant  motion  of  the  surge 
head  is  more  complex  than  that  indicated  by  the  simple  slab 
model  treated  here.  (See  discussion  for  more  details.)  Note  that 
the  slab  velocity  Vtt  can  be  much  larger  the  E  x  B  speed  and  is 
in  a  different  direction  From  this  definition  of  and  equa¬ 
tions  (6)  and  ( 1 1 )  we  find  that 

K.  =  K  {a2R2  +  (1  +  R2(l  -  a))2!1  2  (13) 

l;sing  the  identity  x2  =  ( 1  -  a)2  -  1  +  2a  in  F„,  one  finds  (see 
appendix)  that  the  total  slab  velocity  can  be  expressed  simply 

K.=  K,QIW  +  i/*2!' 2  (14) 

which  is  independent  of  x  and  where  M  is  proportional  to  the 
lotal  h  x  B  speed  as  given  by 

M  =  yd[E02  *  E„2\'2/E0  (15) 

niig  18)  we  have 

M  =  Ml  +  «2(I  -  *>*’!' 2  (16) 

il  has  been  noted  by  Opyenoorth  el  al.  [  1 983]  that  the 
measured  electric  field  within  and  south  of  the  auroral  forms 
...nnected  with  the  passage  of  a  WTS  was  typically  below  the 
thresh- 'id  i  '  15  m\  ml  of  the  STARE  radar  This  gives  an 
-pper  limit  to  the  drift  velocity  M  of  0.375  km's.  For  K. 

qua'iop  1 1 4 n  tv >  equal  the  measured  surge  velocity  of  2.0-  3.3 
mi  -  QH  must  be  therefore,  on  the  order  of  10.  This  implies 
from  figure  2  that  the  mean  electron  precipitation  energy  was 
n  the  order  of  I  keV.  which  is  consistent  with  observations 
n  ;»’n  i ■/  ,./.  1983].  More  comprehensive  auroral  cam- 
faigiiv  a'e  needed  to  simultaneously  define  the  parameters 
s.\!  r.  •  'Ui  model 

-  lint  is  in  f  in  iron  Ion  Recombination 
In  th's  section  we  include  the  effects  of  the  electron  ion 
"no.  h'nation  rate  on  the  WTS  One  exact  solution  to  (4)  can 
:»  •••''!,. ir.cvl  bv  a"..m-Rg  an  ionospheric  magnetospheric  cou- 
v  sv.cl.  itia  .  sc.  '  to  a  constant.  The  equation  for  IH 
..'  the  noril.ern  oounduty  is  now 

<  Lh  1 1  —  —  1„<  IH  t'z  —  GI„  (17) 


<j  =  nrB  ell  (18) 

',,.1  w  tier.  I  „  is  Uelircvl  in  (6|. 

V.  e  now  It.ok  fi  t  a  volution  to  (I7)  that  is  stationary  in  a 
'  '■  tin.oe  vi'icii  i  moving  with  a  velocity  V  relative  to  the 

v.irtlt'.  snrtace  Accordingly,  we  make  a  Galilean  transforma- 
•  "ti  ■  -  I  .  i  I  his  allows  1 1  7]  to  he  rewritten  as 

2.„;  -  G  (I  M  (19) 

h.  .  1,,  2_hi  i  .imt  L„  iT„  C :  Note  by  definition  that 
i  ,9i  ,s  solely  a  function  of  in  the  moving  frame  of  reference 
.  ..o  therefore.  IHi.-  I  is  stationary  in  this  frame  The  solution 

v-h' ~  i  2iM„(t  MH  ~  1*  -  I  (201 

e  2_„.  iHiO|  Hence,  choosing  x  =  constant  implies  that 
’.  Tictiv  tv  profile  is  hyperbolic  in  the  moving  frame.  We 
•  <'  to  fe  where  the  boundary  region  joins  the  main 

::  reg  oi  ,  i  i)  lihester  Baumjohann  model.  Therefore, 


has  the  same  value  as  the  conductivity  inside  the  surge 
head. 

The  assumed  constant  velocity  V  can  now  be  expressed  in 
terms  of  the  conductivity  gradient  at  =  0.  To  see  this  we 
first  take  the  derivative  of  (20)  and  evaluate  it  at  z  —  0  The 
resulting  expression  is 

2-ho  =  ((  -  („)  (21) 

where  £Ho  is  the  value  of  the  conductivity  gradient  at  r  =  0. 
The  physical  meaning  of  (21)  is  made  clearer  by  using  equa¬ 
tion  (2)  and  solving  for  V  We  hnd 

F  =  („  -  GXHo  2 /■.,,*  ;  (0)  (22) 

Wc  have  previously  shown  that  l„  increases  dramatically  as 
the  incident  energy  spectrum  hardens  The  recombination 
term  is  independent  of  Q  and,  therefore.  1  also  increases  in  the 
same  manner. 

Evidently,  the  slab  boundary  will  not  move  poleward  unless 

2,  (0)  >  GT-H„2E0x  f.  =  a,n2e  Q  (23) 

Taking  ar  *  10  '  m  '  s  [  Walls  and  Dunn.  1974],  n-  3 
x  I0M  ions  m\  and  Q  =  5  x  10  4  (eleetron-m)  1  (E,„  =  1 
keV)  we  find  a  threshold  current  of  2.9  jiA  m2  This  value  is 
consistent  with  typical  auroral  current  densities 

The  precipitation  current,  therefore,  must  exceed  some  mini¬ 
mum  value  for  the  surge  to  propagate.  Equation  (22)  also 
implies  that  if  the  current  is  below  the  threshold  value  the 
surge  will  propagate  in  a  negative  or  equatorward  direction 
This,  we  believe,  describes  the  ionospheric  response  during  the 
recovery  phase  of  a  magnetic  substorm.  The  diminishing  pre¬ 
cipitation  current  cannot  sustain  the  high  conductivity  against 
electron  ion  recombination  and  the  surge  head  retreats 

Equation  (17)  can  be  solved  at  the  western  slab  boundary  in 
exactly  the  same  way.  except  now  l„  is  replaced  by  I.  as 
defined  in  (I  I)  The  resul’mg  effective  velocity  depends  on  the 
conductivity  gradient  at  the  western  boundary.  The  conduc¬ 
tivity  gradients,  therefore,  must  also  influence  the  surge's  di¬ 
rection  of  propagation  as  derived  in  section  2 

It  should  be  emphasized  that  1 1 7)  could  be  solved  because 
we  assumed  a  specific  ionospheric  magnetospheric  coupling 
such  that  a  stationary  solution  was  possible  in  a  moving  frame 
of  reference  No  doubt  the  dynamical  nature  of  the  coupling  is 
such  that  the  conductivity  profile  is  probably  nonstationarv  in 
all  frames  of  reference  However,  the  present  theory  qualita¬ 
tively  explains  much  of  the  surge  phenomena  by  using  pres¬ 
ently  available  data 

4  Dis<  t  sshin 

It  has  been  observed  that  the  W  I  S  moves  in  a  series  of 
discrete  steps  or  jumps  [  H  icm  and  Rasioker.  1975.  Pvtie  el 
al .  1976]  This  resull  led  Hininker  er  al  [1980]  to  define  a 
magnetospheric  suhstorm  as  allowing  a  series  of  multiple 
surges  during  the  expansion  phase  Each  surge  corresponds  lo 
an  individual  suhstorm  onset  that  on  a  time  average  shifts  the 
maximum  poleward  expansion  northwestward  In  the  present 
model  each  surge  corresponds  lo  a  temporal  hatdcmng  of  the 
precipitating  electron  flux  energy  spectrum  I  his  temporal 
hardening  is  related  lo  the  suhstorm  onset  processes  in  the 
magnetotail  and  the  generation  of  parallel  electric  fields  which 
is  outside  the  scope  of  the  present  work 

Samson  amt  Rostoker  [1983]  noted  Ihe  presence  of  associ¬ 
ated  Pi  2  signatures  with  the  W  I  S  The  surge  marks  the  Iran 
sition  from  equatorward  to  poleward  type  Pi  2  polarization 
signatures  The  elliptical  polarization  of  Ihe  Pi  2  s  is  con 
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siilered  to  he  caused  by  (he  longitudinal  expansion  of  field- 
al i lined  currents  in  the  surge  together  with  the  conjugate  re¬ 
flection  of  field  aligned  current  pulses.  We  suggest  that  the 
same  enhanced  held -aligned  current  produces  a  jump  in  the 
surge  velocity  as  predicted  by  the  model  presented  here  and  as 
observed  by  U  icns  mid  Raatoker  [  1 9 7 5  J  and  Pvtte  el  al. 
(  1 9  7  (s  | 

The  surge  head  ma>  at  times  break  oil  from  the  main  hodv 
forming  a  separate  surge  (See  Figure  11  in  the  paper  bv 
.Sum si >n  mid  Rostoker  [1983]  )  In  terms  of  the  present  theory, 
once  a  new  conductivity  gradient  is  established  between  the 
surge  regions  their  relative  velocity  is  determined  by  the 
energy  spectra  of  the  precipitating  electrons  at  the  westward 
head  of  the  separate  parts  The  propagation  of  each  part  can 
be  described  by  modeling  them  as  a  separate  slab  Instanta¬ 
neous  velocities  of  up  to  3(1  km  s  have  been  observed  at  the 
leading  branch  of  the  western  electrojet  [Opyennorth  et  ui, 
1983.  (i  Rostoker.  private  communication.  1983]  This 
branch  can  be  considered  as  a  separate  slab  with  Us  own 
model  parameters.  Westward  velocities  of  30  km  s  will  result 
from  (1 1)  if  we  assume  En  =10  2  V  m,  x  =  0.  R  =  3,  and  the 
incident  electron  energy  to  be  approximately  5  keV.  The 
higher  westward  velocities  at  (he  same  precipitation  energy 
are  due  to  the  additional  R  factor  in  l  „,  which  arises  from  the 
enhanced  Cowling  current  when  i  -  0  Yahnin  et  til  [1983] 
note  that  a  harder  precipitation  energy  exists  inside  the  surge 
in  comparison  with  outside  the  surge  region  which  is  consis¬ 
tent  with  this  result 

Surges  have  also  been  occasionally  observed  moving  in  a 
northeasterly  direction  [/’vfte  et  al.,  1976].  This  can  be  ex¬ 
plained  in  the  present  model  by  overclosurc  of  the  Hall  cur¬ 
rent  on  the  poleward  surge  boundary  as  shown  in  Figure  4. 
Intense  electron  precipitation  (a  >  1)  causes  a  poleward  polar¬ 
isation  electric  field  with  an  associated  eastward  Hall  current 
that  weakens  the  westward  Pedersen  current  driven  by  the 
external  electric  held  E„. 

The  direction  of  the  surge  motion  is  related  to  the  direction 
of  the  resultant  electric  field  from  Figure  4  it  is  noted  that  a 
northwestward  direction  corresponds  to  i  i  (I  8  0  9  Using 
this  result  in  (81  for  l.p  implies  that  the  resultant  electric  field 
points  in  the  southwestern  direction  that  is  consistent  with  the 
results  of  Inhester  et  til  [  1981  ] 

Out  icsults  in  section  2  (sommnri/cd  in  I  able  I  lot  x  -  ft 
indicated  that  an  unphysical  sutge  velocity  could  oectu  if  the 
incident  energy  spectrum  is  too  hard  However,  in  section  t 
we  found  that  the  surge  velocity  is  decreased  by  the  electron 
ion  recombination  rate  and  that  the  magnitude  of  the  precipi¬ 
tating  current  is  important  Note  also  that  the  relation  be¬ 
tween  x  and  the  incident  energy  spectrum  is  unknown  so  that 
the  excessive  velocities  given  in  Table  I  probably  do  not 
occur 

In  summary,  it  is  found  that  the  present  model  predicts  a 
wide  range  of  complex  surge  phenomena  depending  on  the 
ratio  of  the  Hall  to  Pedersen  conductivities,  the  degree  of 
ionospheric  current  closure  into  (lie  magnetosphere  and  the 
energy  spectra  of  precipitating  electrons  The  following  sum¬ 
marizes  our  main  results 

1  We  find  that  the  direction  of  the  W  I  S  depends  stronglv 
on  the  degree  of  current  closure  on  the  poleward  boundary 
The  sensitivity  of  the  surge  direction  to  closure  depends  on  the 
ratio  (R|  of  the  Hall  to  (he  Pedersen  conductivities 

2  The  magnitude  of  the  surge  velocity  is  sensitive  to  the 
energy  spectrum  of  the  precipitating  electrons  and  weakly  de¬ 
pendent  on  R  The  ratio  of  the  surge  velocity  to  the  measured 


drift  velocity  is  independent  of  the  degree  of  closure  al  both 
the  surge  head  and  the  northern  boundary 

3  The  expansion  phase  of  the  substorm  is  explained  by 
assuming  that  the  initial  arc  brightening  arises  from  a  sudden 
hardening  of  the  precipitating  electron  energy  spectrum  at  its 
poleward  boundary. 

4  Inclusion  of  electron  ion  recombination  effects  high¬ 
lights  the  role  of  the  precipitating  current  intensity  in  modu¬ 
lating  the  surge  propagation  and  explains  the  equatorward 
retreat  of  the  surge  during  the  substorm  recovery  phase 

5  Details  of  the  surge  propagation  depend  on  how  the 
magnetosphere  and  ionosphere  are  coupled  as  reflected  in  the 
functional  form  of  x  and  Q.  Hence  the  energy  source  is  clearly 
Icoated  in  the  magnetosphere  and  a  complete  description  of 
substorm  phenomena  must  take  this  into  account 

APPI  NDIX 

Fquation  ( 14)  is  obtained  from  ( 1 3)  by  using  the  given  ident¬ 
ity  The  intermediate  results  for  the  bracketed  term  in  (13)  are 

(I  -  a)’/?2  +  I  +  R1  +  R4(l  -  3)1 

=  (I  -  a):R2[I  +  R2]  +  1  +  R2 

=  II  +  R;)[  1  +  (1  -  x) 2R2]  (Al) 

Comparison  of  this  expression  with  (16)  immediately  gives 
(14) 
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A  model  is  developed  lh.it  relates  the  poleward  leaps  of  the  westward  traveling  surge  (WISl  and  the 
generation  of  Pi  2  pulsations  The  feedback  instability  developed  by  Sato  and  coworkers  is  combined 
with  the  dynamic  surge  model  of  Roihwell  et  al  (I9K4I  We  lind  that  our  previous  results  on  the  motion 
of  the  W  I  S  are  related  to  the  zero-order  terms  in  the  Sato  formulation  The  linearized  first-order  terms 
give  rise  to  a  dispersion  relation  with  solutions  in  the  Pi  2  frequency  range  The  /th  frequency,  /,.  is  found 
to  be  turned  on  in  integer  multiples  of  one-half  the  Alfven  bounce  frequency  between  the  ionosphere  and 
the  plasma  sheet  Most  important,  however,  is  that  the  /th  frequency  is  turned  on  only  when  the 
precipitating  electron  energy  exceeds  a  certain  value.  As  previously  shown,  the  velocity  of  the  poleward 
surge  boundary  also  increases  when  the  energy  of  the  precipitating  electrons  is  enhanced  Therefore  the 
poleward  leap  of  the  surge  during  substorm  onsets  is  accompanied  by  the  generation  of  higher  Pi  2 
frequency  components.  The  time  evolution  of  the  composite  Pi  2  pulse  is  obtained  using  the  calculated 
decay  rates,  and  agreement  with  data  is  shown 


In  rRomtcTtoN 

A  recent  review  of  Pi  2  pulsations  and  substotm  onsets  has 
been  given  by  Baumjohann  and  GU tssmeier  [1984],  On  the 
basts  of  previous  work  by  Rostoker  anil  Samson  [1981]  and 
Samson  and  Rostoker  [198.1]  there  is  clearly  an  observed  re¬ 
lationship  between  the  westward  traveling  surge  (WTSl  and 
the  generation  of  Pi  2  pulsations  during  substorm  onsets.  In 
particular,  the  WTS  marks  the  longitudinal  transition  from 
the  equatorward  to  poleward  Pi  2  A  poleward  Pi  2  exists 
within  the  surge  head  and  to  the  east  An  equatorward  Pi  2 
predominates  equatorward  and  to  the  west  of  the  surge  [Ros- 
loker  and  Samson.  1981]  The  maximum  intensities  of  the  Pi  2 
pulsations  were  found  along  the  equatorward  boundaries  of 
the  elcctrojcts  This  led  Rostoker  and  Samson  [1981]  and 
Samson  and  Rostoker  [1981]  to  suggest  that  the  resonance 
region  of  the  Pi  2  pulsations  is  localized  within  the  surge 
region  and  is  consti aitietl  to  lem.im  on  closed  tield  lines  Ros- 
toker  and  Samson  |I9KI|  also  suggest  that  the  llarang  dis¬ 
continuity  is  the  energy  source  region  for  the  Pi  2  Samson  and 
Harrold  [1981],  using  the  University  of  Alberta  magnetometer 
chain,  found  that  within  the  WTS  the  Pi  2  polarization  pat¬ 
terns  arc  clockwise  (CW)  as  viewed  downward  On  the  other 
hand,  lister  el  id.  [1984]  found  with  the  mid-latitude  Air 
Force  Geophysics  Laboratory  (AI  G1.)  chain,  at  the  same  lon¬ 
gitude  but  equatorward.  that  the  polarization  cllipltctty  is  pre¬ 
dominantly  counterclockwise  ((.'(.)  These  results  arc  consis¬ 
tent  with  the  fact  that  equatorward  and  poleward  of  the  WTS 
the  polarizations  arc  counterclockwise  while  far  to  (he  east 
and  west  of  the  WTS  they  are  CW  The  Pi  2  polarizations 
relative  to  the  WTS  are  therefore  quite  complex 
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Causes  of  the  observed  east-west  north-south  transitions  in 
the  Pi  2  polarizations  are  an  open  question.  According  to 
Pashin  et  al.  [1982],  the  Pi  2  current  system  is  primarily  lo¬ 
cated  in  the  upward  current  at  the  surge  head,  and  the  west¬ 
ward  movement  of  the  surge  causes  linearly  polarized  Alfven 
waves  to  appear  elliptically  polarized  on  the  ground.  Samson 
[1982],  on  the  other  hand,  explains  the  observed  polarization 
patterns  in  terms  of  a  longitudinal  distribution  of  oscillating 
current  sheets  Lester  et  al.  [1983.  1984]  interpret  the  mid- 
latitude  Pi  2  pulsations  as  arising  from  the  field-aligned  cur¬ 
rents  that  form  the  substorm  current  wedge  Ellis  and  South- 
wood  [1983]  have  examined  the  reflection  of  Alfven  waves 
from  ionospheres  with  a  discontinuity  in  either  the  Hall  or  the 
Pedersen  conductivity  The  reflection  properties  of  the  Alfven 
wave  depend  not  only  on  the  discontinuity  type  but  also  on 
the  orientation  of  the  incident  electric  field  vector  In  two  of 
the  four  eases  studied,  the  field-aligned  current  sheets  act  as 
subsidiary  surface  waves  centered  on  the  field  lines  connected 
to  the  discontinuity.  The  subsidiary  waves  are  circularly  polar¬ 
ized  and  suppress  any  net  flow  of  Hall  current  across  the 
discontinuity  Glassmeier  [1984]  has  extended  the  work  of 
Tilts  and  Southwood  to  include  arbitrary  distributions  in  the 
height-integrated  conductivity.  Recently.  Southwood  and 
1 1  lathes  [1985]  have  suggested  that  two  oppositely  traveling 
easi-west  surface  waves  parallel  to  the  conductivity  gradients 
cotilv!  give  a  combined  signal  that  reproduces  many  of  the 
observed  features  of  the  Pi  2  pulsations  An  important  unre¬ 
solved  question  is  the  relation  of  the  Pi  2  current  system  to  the 
overall  current  system  that  forms  the  substorm  current  wedge. 

It  ts  assumed  that  the  initial  diversion  of  the  cross-tail  cur¬ 
rent  is  carried  by  a  transverse  Alfven  wave  f Baumjohann  and 
Glassmeier.  19X4]  The  subsequent  reflection  of  the  wave  [Ma- 
Ifmekrodt  and  (orison.  |9"*S.  \rshtdo.  1979.  Kan  el  al..  1982] 
between  the  conjugate  ionospheres  and  the  triggering  of  sec¬ 
ondary  Allven  waves  in  the  ionosphere  [  Maltsev  et  al..  1974; 
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1  ig  i  Idealized  model  of  llie  W  I  S  as  proposed  by  Inhesier  i'I  al  |  19X1  j  1  he  lol.d  external  eleetne  Irehl  /  „  drives  a 
westward  Pedersen  and  a  northward  It. til  current  If  the  northward  Hall  current  is  not  fully  vonftmted  In  licld-.digried 
eurrenls  into  the  magnetosphere,  polarization  charges  hutltl  up  on  ihe  poleward  surge  boundary,  producing  a  southward 
directed  electric  field  £p.  This  southward  electric  field  produces  a  Hall  current  in  the  same  direction  as  Ihe  Pedersen 
current  from  ihe  original  electric  held 


Taman  anil  Mmra.  1982]  both  contribute  to  the  Pi  2  pulsation 
tune  profile  As  pointed  out  by  Baionjohann  and  Glassmeier 
[19X4],  it  is  important  to  treat  (he  ionosphere  as  an  active 
part  of  the  tonosphere-magnetospheric  coupling  rather  than 
just  a  reflecting  boundary 

I  he  purpose  of  this  paper  is  to  relate  Pi  2  pulsations  to  the 
dynamics  of  the  westward  traveling  surge  (WTS).  A  dynamical 
surge  model  has  been  developed  by  Rolhwell  el  at.  [1984] 
using  the  Inhester-Baumjohann  [ Inliesler  el  at ,  1981 ;  Baumjo- 
ham i,  1983]  representation.  Here  the  feedback  instability 
analysis  [Sum  and  llolzer,  1973;  Holzer  and  Sam ,  1973;  Sam, 
1982]  is  applied  to  the  Inhester-Baumjohann  model  of  the 
WTS  used  by  Rolhwell  el  al.  [1984],  hereinafter  called  paper 
1  The  equations  are  linearized  in  the  standard  fashion,  and  it 
is  shown  that  the  zero-order  terms  recover  the  poleward 
motion  of  the  surge  boundary  as  derived  in  paper  1.  The 
first-order  terms  give  rise  to  a  dispersion  relation  for  the  feed¬ 
back  instability.  This  relation  is  solved  for  the  allowed  fre¬ 
quencies  and  their  associated  growth  rates.  It  is  found  that  the 
number  of  frequencies  generated  is  related  to  the  speed  of  the 
WTS  boundaries.  The  resulting  composite  pulse  shapes  arc 
shown  to  be  very  similar  to  those  measured  by  Singer  et  al. 
( 1985]  at  mid-latitudes.  The  physical  picture  presented  is  that 
the  initial  precipitation  caused  by  the  onset  in  the  magnetotail 
triggers  the  feedback  instability  in  the  coupled  ionosphere - 
magnctosphcriv  s;  ‘  'em,  thereby  producing  Pi  2  pulsations 
[ Baumjohunn  and  Glassmeier,  1984].  In  the  next  section  we 
relate  the  feedback  instability  to  the  Inhester-Baumjohann 
WTS  model 

Tut.  Inhlsti  r  Baumjoiiann  Model  and 
t hi.  Feedback  Insiaweiiy 

The  Inhester-Baumjohann  model  is  shown  in  Figure  I.  A 
uniform  westward  electric  field  F.0  produces  a  northward  Hall 
current  that  closes  into  the  magnetosphere  via  field-aligned 
currents  along  the  poleward  surge  boundary.  If  the  precipi¬ 
tation  is  insufficient  to  close  oil'  this  Hall  current,  then  an 
eflective  polarization  charge  builds  up  along  the  conductivity 
gradient,  producing  a  southward  directed  electric  field  This 
electric  field  creates  a  southward  Pedersen  current  opposite  to 
the  original  Hall  current  and  a  westward  Hall  current  that 
adds  lo  the  original  westward  Pedersen  current  driven  by  li„. 

As  noted  above,  the  Pi  2  current  system  and  the  substorm 
current  system  inside  the  wedge  are  not  always  the  same  In 
this  paper  we  initially  assume  that  both  current  systems  are 
colocated,  which  is  true  approximately  65%  of  the  time 
(  /  osier  el  al .  198  '  | 

Ihe  feidback  instability  [ (h/awa  and  Sam.  1971;  Sam  and 
lie. ri.  IV  3  his  iis  foundation  in  earlier  work  bv  Atkinson 


[1970],  The  feedback  instability  was  originally  developed  for 
the  quiet  arc,  which  has  a  much  larger  extension  in  the  cast 
west  direction  than  in  the  north-south  direction,  similar  to  the 
surge  region.  A  downward  field-aligned  current  on  the  cqu.i- 
torward  edge  of  the  arc  closes  via  a  northward  Pedersen  cur¬ 
rent  lo  an  upward  field-aligned  current  on  (he  polcwaid 
boundary.  According  to  Atkinson  [1970]  a  local  ionospheric 
conductivity  enhancement  causes  a  local  decrease  of  the  elec¬ 
tric  field.  The  resulting  divergence  of  the  maenctospheric  po¬ 
larization  currents  produces  precipitation  that  increases  the 
original  conductivity  enhancement.  If  one  visualizes  a  north- 
south  ionospheric  wave  in  the  model  in  Figure  I.  then  one  will 
have  periodic  conductivity  enhancements  which  could  lead  to 
multiple  arcs.  The  analysis  for  the  quiet  arcs  as  given  by  Sam 
and  Holzer  [1973]  and  Holzer  and  Sam  [1973]  is  similar  to 
the  analysis  presented  here  In  their  analysis,  active  and  pas¬ 
sive  ionospheric  regions  are  conjugately  connected  by  the 
same  field  line.  The  active  ionosphere  acts  as  an  ac  generator 
which  produces  an  Alfven  wave  that  is  damped  in  the  passive 
ionosphere.  Sum  [1978.  1982]  dispensed  with  the  conjugately 
connected  active  and  passive  ionospheres  and  required  that 
the  Alfven  waves  reflect  upon  reaching  the  equatorial  plane 
Sam  |  1978]  also  notes  that  the  theory  of  quiet  arcs  must  be 
self-consistent  with  the  presence  of  a  westward  electric  field. 
This  is  also  a  feature  of  the  WTS  as  seen  in  Figure  I 

The  feedback  instability  works  because  electrons  tend  to 
(low  toward  the  positive  part  of  the  potential  perturbation 
along  the  field  lines  The  inductive  reactance  of  the  mag¬ 
netosphere.  however,  causes  a  phase  lag  in  the  precipitation 
such  that  it  adds  to  the  original  ionization  enhancement,  caus¬ 
ing  the  instability  to  grow  Should  the  magnetosphere  have  a 
capacitativc  reactance,  then  the  precipitating  electron  flux  co¬ 
incides  with  the  valley  of  the  density  distribution,  and  the 
perturbation  decays.  The  magnetosphere  must  have  an  induc¬ 
tive  response  for  the  feedback  instability  to  occur.  By  induc¬ 
tive  and  capacitativc  reactance  we  refer  to  the  eflective  termi 
naled  transmission  line  impedance  for  Alfven  waves  along  the 
magnetic  field  lines  [ Sam.  1982]  T  his  work  has  been  extended 
by  Mmra  and  Sam  [1980]  lo  the  global  formation  of  multiple 
auroral  arcs. 

Taman  and  Mnira  [1982],  Malta  el  al  [1982].  and  Taman 
[1984]  considered  a  nonunilorm  magnetosphere  coupling  to 
the  ionosphere  Negative  Joule  dissipation  in  the  ionosphere  is 
accompanied  by  a  growing  oscillation  and  an  outflowing 
Poynling  (lux  from  the  ionosphere  Damped  oscillations  occur 
when  energy  is  supplied  from  the  magnetosphere  to  Ihe  iono¬ 
sphere  A  large-scale  umfoim  elect nc  field  in  Ihe  ionosphere  is 
required  to  drive  the  grow  mg  instability  ' 

In  the  present  work  the  generation  ol  Pi  2  pulsations  is 
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considered  as  follows.  The  injection  of  hot  electrons  from  the 
plasma  sheet  during  a  substorm  onset  initiates  the  feedback 
instability  in  the  ionosphere-magnetosphere  system.  This  as¬ 
sumption  is  based  on  the  satellite  observations  of  Sakurai  and 
XtcPIwron  [1983]  that  the  Pi  2  burst  is  superimposed  on  a  dc 
shift  in  the  azimuthal  component  of  the  magnetic  field  which 
is  caused  by  field-aligned  currents  From  their  Figure  18  the 
Pi  2  occurs  at  the  beginning  of  the  dc  shift  The  parallel  (pre¬ 
cipitating)  current  is  therefore  considered  to  be  decomposable 
into  ac  and  dc  components  The  dc  component  is  the  primary 
injection  front  the  distant  magnetotail,  and  the  ac  component 
arises  from  the  initial  injection  transient  and  the  feedback 
instability.  It  will  be  shown  below  that  the  number  of  modes 
that  are  stimulated  is  dependent  on  the  energy  of  the  zero- 
order  precipitation  It  is  argued  that  the  feedback  instability 
does  not  reach  the  nonlinear  stage  examined  by  Sato  [1978] 
and  Miura  et  al.  [1982]  since  the  flux  associated  with  ;  (dc) 
raises  the  ionization  level  within  a  time  u  (o,.Y„)  where  n, 
is  the  electron-ion  recombination  rate  and  V0  is  the  zero- 
order  ion  density.  After  this  time,  electron-ion  recombination 
dominates,  and  the  composite  pulsation  decays 

Kan  t’t  al.  [1982]  consider  the  Pi  2  wave  form  to  arise  from 
the  superposition  of  the  reflected  and  incident  Alfven  waves 
impinging  on  a  passive  ionosphere.  T  hey  neglect  polarization 
and  Hall  current  effects.  Here,  on  the  other  hand,  we  consider 
the  natural  modes  arising  from  a  self-consistent  ionosphere- 
magnetosphere  interaction  (i.e..  the  feedback  instability). 

Formula  nos  or  tut  Model 

In  paper  1  [Rothwell  et  al.,  1984]  the  solutions  to  the  time- 
dependent  zero-order  equations  along  the  surge  boundary 
conductivity  gradients  were  solved  On  the  other  hand,  in 
Sato's  [1982]  theory  the  zero-order  enhanced  ionization  den¬ 
sity  Na  is  taken  as  constant  in  space  and  time,  which  is  consis¬ 
tent  with  the  Inhcster-Baumjohann  model  inside  the  surge 
region. 

In  formulating  the  present  Pi  2  pulsation  model  we  first 
connect  the  magnetospheric  transverse  Alfven  wave  with  the 
ionospheric  drift  wave  [7  amao  and  Miura,  1982].  The  Alfven 
wave  is  probably  kinetic.  However,  the  perpendicular  scale 
lengths  considered  here  are  much  larger  than  the  ion  gy- 
roradius  so  that  the  Alfven  dispersion  relation  is  essentially 
the  same  as  in  the  MHD  case.  It  is  permissible,  therefore,  in 
the  present  context  to  use  the  MHD  dispersion  relation.  An 
Alfven  wave  carries  parallel  current  which  is  related  to  the 
divergence  of  the  transverse  electric  wave  field  by 

7.„j  -  V  •  £  (1) 

where  Ztl  is  the  characteristic  impedance  of  an  equivalent 
transmission  line  terminated  at  both  ends  by  the  impedance 
/„  and  /  is  in  the  same  direction  as  the  ambient  magnetic 
field  It  can  be  shown  [see  Fan  et  ill ,  1982]  that  j  and 
V  ■  £  satisfy  the  transmission  line  equations  Sato  [1982] 
used  the  transmission  line  analogy  to  impose  the  ionospheric 
and  equatorial  boundary  conditions  As  seen  at  the  iono¬ 
sphere,  the  magnetospheric  impedance  is  given  by 


The  northward  ionospheric  current  component  in  the 
Inhester-Baumjohann  model  is  given  by 

J,  =  £,!„  +  £,  I,  (3) 

where  £„  and  are  the  Hall  an  Pedersen  conductivities, 
respectively,  and  £  is  the  electric  field.  The  conductivities  are 
normalized  to  their  zero-order  values  by 

(4) 

S,  = 

where  and  lpl,  are  uniform  inside  the  surge  region  and  .V 
is  the  height-integrated  ionospheric  ionization  density  normal¬ 
ized  to  the  uniform  density  within  the  surge  region  [Sato, 
1982]  Now  we  similarly  have  for  the  westward  (y|  current 
component 

J .  =  v,  -  (51 

Equations  (3),  (4).  and  (5)  are  linearized  as  follows. 

E.  =  F.,0  +  Ex 

£,.  =  £,.  o  +  £y  (6) 

A-  =  I  +  N 

where  the  values  with  overbars  are  of  first  order.  £>0  is  the 
primary  electric  field  that  drives  the  substorm  current  wedge, 
and  £>f)  is  the  primary  north-south  polarization  electric  field 
Both  £y0  and  £l0  are  assumed  constant  inside  the  surge 
region.  Now  the  constancy  of  N„  and  £0  inside  the  surge 
region  implies  that  there  is  no  net  zero-order  field-aligned 
current  closing  into  the  magnetosphere  inside  the  surge.  How¬ 
ever,  there  is  still  a  zero-order  energetic  electron  precipitation 
that  maintains  the  high  conductivity  inside  the  surge  region. 
The  current  carried  by  the  energetic  electrons  must  be  precise¬ 
ly  balanced  bv  a  flux  of  upward  flowing  ionospheric  electrons 
of  lower  energy.  In  this  manner  a  high  conductivity  level  is 
maintained  with  no  net  current  closure.  The  divergence  of  J  is 
given  by 

V  •  3  =  J„  •  VN  +  Ir0Y  ■  E  +  Imi(Y  x  E),  (7) 

where  J0  is  the  zero-order  two-dimensional  ionospheric  cur¬ 
rent.  The  ionosphere  is  considered  a  source  of  transverse 
Alfven  waves  so  that  (V  x  £).  =  0  [Malice v  et  a!.,  1974]  ex¬ 
cludes  magnetoacoustic  waves  in  the  ionosphere.  The  first- 
order  field-aligned  current  density  is  given  by 

I  -  V  •  J  -  -  J0  •  YiV  -  Ip0V  •  E  (8) 

We  assume  that  the  first-order  current  is  being  carried  by  hot 
electrons  Note  that  we  have  a  coordinate  system  with  the 
upward  current  as  positive,  which  means  that  we  have  a  minus 
sign  in  the  following  relation. 

/,  -  -v  •  f:  (9i 

and  upon  combining  (8)  and  (9)  wc  have 

j,  =  -Jo-VNll  -ZIp0>  (10) 


/  //,,  col  tod  l \)  (2) 

where  f is  the  magnetospheric  Alfven  speed  and  I  is  the 
length  of  the  field  line  between  the  ionosphere  and  the  equa¬ 
tor  This  expression  also  assumes  that  ;  is  zero  at  the  equa¬ 
torial  plane  7.,,  is  given  by  n„VA.  Then  (II  holds  at  the  iono¬ 
sphere  if  /  in  ( 2)  replaces  in  ( 1 1. 


Now  the  first-order  continuity  equation  gives  us 
fiiY  el  -  -Qhj  eSa  -  2 i,.V 
5  =  a,  \  o  h 


(111 


where  Qh  is  the  height-integrated  ion  production  efficiency 
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[Rees,  1963],  The  insertion  of  ( 10)  into  (II)  gives 

r.V  rt  =  -  V  •  V ,*?'(!  -  ZI^)  -  2*A?  (12) 

where  the  components  of  V  are  given  by 

.  (13) 

l,  =  (HiFJl  +  RJ(I  -  jt)],R 

The  parameter  x  is  the  closure  parameter  as  defined  in  the 
paper  by  Rothwell  et  al.  [1984],  is  the  1x8  drift  velocity 
taken  as  *0.25  km  s,  and  R  is  the  ratio  of  the  height- 
integrated  Hall  to  Pedersen  conductivities,  which  we  take  to 
be  *  3.  These  are  the  components  of  the  surge  velocity  as 
derived  by  Rothwell  et  al.  [1984]  without  the  electron-ion 
recombination  term.  This  velocity  term  was  derived  in  (12)  by 
noting  that 

QhJ0  eS0  =  Qh  e/V„(E  •  E)  (14) 

where  X  *  vN„  B„  leads  to  equations  (9)  and  (1 1)  in  the  paper 
by  Rothwell  el  al.  [1984], 

It  is  assumed  in  the  usual  manner  that 

V  *  exp  [i(k,x  +  kty  -  to/)]  (15) 

where  <■>  is  considered  complex  (w  =  to,  +  ito, )  and  k  is  the 
wave  vector  for  the  ionospheric  wave.  Taking  X  =  'S.pUZn  cot 
[i  irl  f ,  -  mt],  we  have 

-no,  +  to,  —  — tV  •  k(l  +  iX)/(l  +  A" 2)  —  2a  (16) 
which  leads  to 

to,  =  V  •  k  II  +  A'2)  (17) 

foi  the  frequency  dispersion  and 

to,  =  to,X  -  2<r,Nn  (18) 

for  the  growth  rate. 

Suit i  [1978,  1982]  has  argued  that  the  maximum  growth 
rate  occurs  for  ,Y  =  I  However,  X  is  not  a  free  parameter, 
and  its  value  must  be  consistent  with  the  solutions  to  (17).  For 
each  value  of  n  In  -  0,  1 ,  2,  ■  •)  one  obtains  two  solutions  to 
ii~'  One  root  cctrevp'nds  to  a  negative  transmission  line 
irpedanee  or  a  capacitative  reactance  and  is  highly  damped 
The  other  root  corresponds  to  a  positive  (inductive)  reactance 
md  gives  positive  growth.  This  result  is  consistent  with  the 
view  ,  of  Sato  [  1978]  In  the  following  discussion  we  consider 
only  the  inductive  roots  of  to,. 

One  v  in  view  the  substorm  current  wedge  as  forming  a 
box  shaped  region  of  enhanced  conductivity.  The  ionospheric 
waves  associated  with  the  Pi  2  pulsations  partially  reflect  off 
the  conductivity  gradients  on  the  boundaries  forming  standing 
waves  m  the  a  and  i  directions.  In  the  following  calculations  it 
is  .  niiied  that  the  waves  arc  dominated  by  the  fundamental 
modes.  ie.  the  wavelengths  in  each  component  are  of  the 
same  order  as  the  scale  size  in  that  direction 

\ow  the  k  ■  V  term  in  ( 1 7)  can  be  expressed  as  (see  ( 1 3)) 

d  z.,k  •  V  \lnQin<)  -  !/i  +  [I  f  K2(l  -  *)]  R\  (19) 

where  -■  is  the  length  of  the  surge  wedge  in  the 

west  direction,  and  /,  is  the  surge  width  in  the  north- 
..  it1  diction.  Figure  2  shows  a  plot  of  the  right-hand  side  of 
’  '  i  various  values  of  a  and  R.  The  long  dashes  in  this 

,  ,  r  c.  in  the  product  yx.  Note  that  for  most  cases,  yx  is 
.liable  approximation,  which  means  that  except  for  very 
.  /eio  clo.'.ri  i  the  north-south  surge  dimensions  domi- 
V  Pi  2  frequency  characteristics  In  the  following  exam- 


Fig  2.  The  parameter  W  as  defined  in  equation  (19)  In  the  text 
This  figure  shows  that  except  for  almost  zero  closure  lx  =  0)  that 
If  -  yx  idenoted  by  solid  lines)  is  a  reasonable  approximation.  This 
graph  indicates  that  the  north-south  surge  dimensions  generally  de¬ 
termine  the  Pi  2  frequency  characteristics 

pie,  therefore,  we  set  /.y  to  infinity  and  treat  the  north-south 
case.  Note  also  that  for  larger  x  the  Pi  2  pulsation  frequencies 
are  higher.  The  dependence  of  W  and  R  is  weak. 

As  noted  above,  we  solve  for  the  inductive  (.Y  >  0)  roots  of 
(17).  The  resulting  values  of  to,  are  inserted  into  the  first  term 
on  the  right-hand  side  of  (18)  and  plotted  in  Figure  3a.  The 
characteristic  north-south  dimension  of  the  surge  is  taken  as 
500  km,  and  the  quantities  E0Z0  and  VA  I  are  set  to  10  and 
0.005  s  ',  respectively.  Note  that  a  mode  is  excited  whenever 
X  =  0  or  to,  =  In  +  \)nVAil.  This  means  that  (17)  has  only 
physical  solutions  for  frequencies  less  than  itf,  kx.  Now  the 
value  of  Vx  increases  with  the  energy  of  the  precipitating  elec¬ 
trons  [Rothwell  et  al.,  1984].  Thus  the  energy  of  the  zero-order 
component  Idc)  of  the  precipitating  electrons  controls  the 
number  of  excited  modes  as  seen  from  Figure  3a.  More  ener¬ 
getic  precipitation  is  associated  with  higher  magnetic  activity. 
Sakurai  and  McPherron  [1983]  analyzed  Pi  2  magnetic  ac¬ 
tivity  Sakurai  anil  McPherron  [  1983]  analyzed  Pi  2  magnetic 
pulsations  observed  at  geosynchronous  orbit  on  ATS  6  They 
found  that  as  magnetic  activity  increased,  the  frequency  spec¬ 
trum  became  more  complex  with  more  spectral  power  at 
higher  frequencies  This  is  consistent  with  Figure  3d  in  that  a 
larger  ),  is  associated  with  more  energetic  electron  precipi¬ 
tation  and  faster  motion  of  the  surge  boundaries.  As  seen  from 
Figure  3d.  higher-frequency  modes  are  excited  at  higher 
growth  rates,  implying  a  greater  contribution  to  the  spectral 
power.  Figure  3h  shows  the  corresponding  mode  frequencies 
for  the  same  inputs  as  given  for  Figure  3a. 

A  horizontal  line  in  Figure  3a  would  represent  the  damping 
due  to  electron-ion  recombination.  The  net  decay  is  the  differ¬ 
ence  between  this  curve  and  the  individual  to,. Y  curves.  The 
present  theory  is  applicable  when  the  intersection  of  the  l'v 
line  with  the  electron-ion  loss  rate  is  above  the  growth  curves 
If  the  intersection  is  below  the  growth  curves,  then  continuous 
growth  is  predicted,  which,  of  course,  is  unrealistic.  Therefore 
this  latter  case  should  be  treated  by  a  more  sophisticated  non- 
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Fig  36  Frequencies  of  the  various  modes  for  the  same  inpul  pa¬ 
rameters  .is  shown  in  Figure  G 


AFGL  MAGNETOMETER  NETWORK  JUT  6,1979 
HIGH  PASS  FILTERED  00033  Hi 


1  lli.!.,..  ,  .  _ _ _ 

c’  to'  O  10° 

v<n. 

Fig  3d.  Growth  rates  for  ZPZV  and  VA  I  =  0005  The  north- 
south  dimension  of  ihe  surge,  a,,  is  taken  as  500  km  The  growth  rates 
m,X  of  the  various  excited  frequency  modes  are  plotted  as  a  function 
of  k’,  a,  where  kj  is  the  aero-order  poleward  surge  velocity,  neglecting 
electron-ion  recombination  effects  The  value  of  this  ratio  determines 
the  number  of  modes  excited  and  the  overall  time  profile  of  the  re¬ 
sulting  Pi  2  pulsation  A  horizontal  line  equal  to  2 <r,n  would  represent 
wave  damping  due  to  electron-ion  recombination 


linear  approach  The  sensitivity  of  the  results  to  the  various 
parameters  is  as  follows.  Higher  electron  precipitation  energy 
implies  that  higher-order  modes  are  excited  with  faster  growth 
rates  and  also  that  a  higher  ionospheric  ionization  density  is 
attained.  These  two  effects  tend  to  offset  each  other  in  the 
sense  that  a  higher  ionization  level  implies  faster  electron-ion 


Fig  4  AFGL  magnetometer  filtered  Pi  2  data  taken  on  July 
1979  (courtesy  H.  Singer  AFGl-f 


16. 


recombination  and  hence  enhanced  damping.  The  faster 
growth  rates  dominate,  however,  so  that  the  net  effect  is  larger 
Pi  2  pulsations  at  higher  incident  electron  energies  Higher 
incident  flux  at  fixed  incident  energy  can  lead  to  overdamping 
lower  flux  causes  more  rapid  growth.  Smaller  values  of  P,'a.t 
(longer  held  lines  I  lead  to  the  excitation  of  more  modes  shifted 
to  lower  frequencies. 

in  the  present  model  an  external  condition  is  needed  to 
relate  the  zero-order  electron  precipitation  flux  with  the  elec¬ 
tron  precipitation  energy  in  order  to  ensure  damped  Pi  2  pul¬ 
sations  We  therefore  took  the  results  of  Fridman  and  l.emaire 
[19K0],  who  relate  the  field-aligned  electron  fluxes  with  the 
associated  field-aligned  potential  drops.  They  consider  five 
separate  cases  corresponding  to  different  boundary  conditions 
in  the  plasma  sh  et  source.  It  was  found  that  aJ  five  cases 
gave  values  for  electron-ion  recombination  damping  that  were 
above  the  growth  rate  curves  shown  in  Figure  3u.  The  h '■id- 
man  and  l.i  tnairi -  [1980]  results  therefore  are  consistent  with 
damped  Pi  2  pulsations  as  derived  from  the  present  model 

Now  in  equation  (231  of  paper  1  we  found  that  ihc  incident 
flux  along  the  poleward  boundary  had  to  exceed  some  critical 
value  in  order  for  the  surge  to  propagate  It  turns  out  that  this 
critical  flax  exceeds  the  required  flux  level  inside  the  surge 
region  to  cause  damped  Pi  2  pulsations  Therefore  from  the 
present  work  one  expects  damped  Pi  2  pulsations  associated 
with  poleward  surge  movements  This  is  an  important  consist¬ 
ent  test  between  our  theoretical  approach  and  observational 
results 

Figure  4  shows  a  Pi  2  pulsation  as  measured  by  Singer  el  a! 
[lugs]  aru|  Figure  6  shows  the  results  of  the  present  calcula- 
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Fig.  5  f  he  result m>>  model  Pi  2  pulsation  lime  profile  based  on  the  example  shown  in  Figure  hi.  This  is  a  linear 
superposition  of  i he  individual  modes  assuming  all  are  initially  excited  with  an  equal  but  arbitrary  amplitude  Figures 
s.i  SI  are  the  results  for  incident  energies  of  2  keV,  4  keV,  6  keV,  and  ft  keV,  respectively.  The  corresponding  precipitating 
electron  fluxes  used  are  .'.5  «  Ilf.  7  *  10  .  I  I  s  10".  arid  1.5  x  10“  (cm *  si  1  The  equilibrium  ionization  level  N„  is  given 
by  ,V,,  if)/  a, l1  '  where  Q  is  ihe  ion  production  rale  (Rees,  1963]  and  a,  is  the  electron-ion  recombination  rale  The 
damping  term  is  given  by  2rr,.V„  so  lhai  a  minimum  value  of  /  is  required  in  order  to  exceed  Ihe  growth  rales  given  in 


Figure  3u  Lower  values  of  /  lead  to  continued  growth  and 
theory . 

L.in  where  the  ..'.'Hi  l  modes  were  considered  initialized  at 
equal  amplitudes  and  all  having  a  sinusoidal  dependence.  The 
four  cases  I'u  3 Jf  cot  respond  to  incident  energies  of  2  keV.  4 
keV.  6  keV  and  8  keV  The  incident  electron  fluxes  used  in 
Figures  Si  SI  were  3  5  x  Hi  .  7  x  10'.  I  I  x  I ()“,  and  15 
>  10"  lent'  si  :.  respectively  Note  that  the  higher-frequency 
pulsations  occur  at  higher  incident  energies  The  agreement 
with  M‘n,i  r  cl  ill  s  |  id, St]  data  is  seen  to  he  quite  good  for 
t  mure  5.  I  he  other  cases  (particularly  Figures  56  and  5J|  are 
i.oi  1 1 1 c o n s i s i c i ■  i  wiilt  obseivaiions,  although  they  are  a  little 
'  >o  icgiil.it  and  List  too  long.  A  slight  increase  in  ionisation 
.election  precipitation  tluxt  would  significantly  decrease  the 
i’-nsaiiot;  dm. (Hi m  li  is  concluded  that  the  different  damping 
ah  lot  the  ticqiiencies  arising  Irom  the  Sato  formulation  can 
'c. ui  ;■  composite  pulses  which  arc  very  sinulai  to  the  mea¬ 
sured  Pi  2  pulsations  in  shape  and  time  duration  Therefore 
i 1 . i . .  i s j > 1 1 c 1 1 ,  gcncialioii  ol  Allven  waves  could  provide  the  pri- 
n.ir  v  sienaturc  'or  Pi  2  pulsations  The  damping  of  the  Pi  2's 
..  1  erv  sensitive  to  increases  in  the  precipitation  current  that 

...I  ..ICC  S  V 

I  3.  ;  c  -uiting  plivsic.il  model  is  as  follows  The  interruption 
Li  u  dt..k  >.  ncnl  in  the  plasma  sheet  causes  electron 
.  c  a,  l.'.ips.  ol  i.ul  held  lines  to  a  more  dipolar 

c  .  a1  .  a  'i,'  i  mation  <>f  a  suhstoi m  current  wedge 

■iic'h  i In-  i.  i  -phcic  The  impact  of  the  initial  precipitation 


to  the  nonlinear  regime  that  is  not  covered  in  the  present 

which  has  as  a  free  energy  source  the  cast-west  electric  field 
£,,,.  in  the  substorm  current  wedge.  The  feedback  instability 
(ills  the  (lux  lubes  from  Ihe  ionosphere  with  Alfven  waves 
which  form  standing  waves  between  conjugate  ionospheres. 
Since  the  injected  electrons  are  presumably  on  closed  held 
lines,  there  is  a  simultaneous  launching  of  Alfven  waves  into 
the  magnetosphere  from  comparable  locations  in  the  two 
ionospheres.  The  ionospheric  conductivity  is  high  inside  ihe 
surge  so  that  waves  once  injected  into  the  magnetosphere  are 
efficiently  trapped  between  the  conjugate  ionospheres  [llutilws 
and  Southwood.  1976]  Mow  do  these  standing  Alfven  waves 
decay 

Sakurai  and  Mi  Pharr  on  [198.7|  note  Pi  2  polarization  re¬ 
versals  in  space  similar  to  those  observed  by  ground-based 
stations  I  hey  also  point  out  that  Pi  2's  in  space  have  a  large 
comprcssional  component  and  therefore  can  propagate  across 
held  lines  as  fast-mode  hydromagnctic  waves  f.STm/er  <■/  til.. 
1983]  The  propagation  speed  is  faster  than  the  Alfven  speed 
since  the  fast-mode  phase  velocity  is  given  as  [  Akhiczer  cl  til.. 
1975] 

L.  =  |F,:  r  f\:|'  ;  (201 

where  (  ,  is  the  sound  speed  It  is  argued  that  these  waves  still 
have  a  velocity  component  parallel  to  the  magnetic  field  so 
that  they  impinge  on  the  ionosphere  outside  the  surge  region 


ionosphere  triggers  ihe  feedback  instability  Outside  the  surge,  however,  the  ionospheric  conductivity  is 
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Pi 2  POLARIZATION  PATTERN 


Fig.  6  Illustrates  a  possible  explanation  for  Ihe  different  Pi  2  polarization  observed  inside  and  oulside  Ihe  subslorm 
current  wedge.  An  incident  Alfven  wave  from  the  magnetosphere  stimulates  the  feedback  instability,  as  described  in  the 
text.  The  subslorm  wedge  acts  as  an  ionospheric  source  for  Pi  2  pulsations  with  clockwise  polarization  Farther  out  in  Ihe 
magnetosphere  these  couple  to  compressional  waves  that  propagate  perpendicular  to  b  Outside  the  substorm  wedge 
region  these  compressional  waves  couple  to  Alfven  waves  that  are  easily  absorbed  by  the  ionosphere  in  regions  of  low 
conductivity  As  viewed  from  the  ionosphere,  these  incident  waves  will  have  counterclockwise  polarization  consistent  with 
observations. 


much  lower,  and  therefore  wave  reflection  is  much  less  ef¬ 
ficient  [Hughes  and  Southwood,  1976;  Ellis  and  Southwood. 
1983].  The  magnetospheric  Alfven  waves  in  these  regions  are 
rapidly  damped  by  Joule  heating  [Hughes  and  Southwood. 
1976]. 

The  Alfven  wave  created  by  the  Sato  feedback  instability 
originates  in  the  ionosphere  and  propagates  into  the  mag¬ 
netosphere.  This  wave  is  assumed  to  have  clockwise  polariza¬ 
tion  looking  into  the  wave  or  counterclockwise  polarization 
looking  in  the  direction  of  propagation.  It  is  also  assumed  that 
the  sense  of  polarization  is  not  affected  by  cross-field  propaga¬ 
tion  via  the  fast-mode  hydromagnetic  wave.  In  the  low- 
conductivity  region  outside  the  substorm  current  wedge  the 
incident  counterclockwise  wave  impinges  on  the  ionosphere 
from  the  magnetosphere.  If  most  of  the  incident  wave  in  this 
region  is  absorbed  by  the  ionosphere,  then  the  polarization 
looking  down  is  counterclockwise.  In  the  high-conductivity 
region  the  launched  waves  arc  reflected  in  the  conjugate  iono¬ 
sphere  and  return  to  the  source  region  where  they  are  highly 
reflected  Therefore  the  sense  of  polarization  in  the  highly  con¬ 
ductive  region  is  determined  by  the  ionospheric  source  charac¬ 
teristic  of  the  feedback  mechanism  rather  than  by  waves  im¬ 
pinging  from  the  magnetosphere.  These  concepts  are  illus¬ 
trated  in  Figure  6. 

In  summary,  the  sense  of  polarization  for  ionospheric 
sources  and  sinks  of  Alfven  waves  looking  down  al  the  iono¬ 
sphere  should  be  reversed  This  interpretation  is  consistent 
with  the  observations  of  Samson  and  Harrnld  [1983]  as  de¬ 
scribed  in  the  introduction 

The  model  has  the  following  features 

I  The  Pi  2  burst  is  a  result  of  the  sudden  diversion  of  the 
tail  current  to  the  ionosphere  [Sakurai  and  \U  Pherron.  1983] 

2.  The  feedback  instability  is  an  ionospheric  source  for 
Alfven  waves.  Larger-amplitude  components  are  generated  at 
higher  frequencies  for  more  energetic  precipitation 

3.  Standing  Alfven  waves  are  created  between  conjugate 
ionospheres  on  field  lines  that  connect  the  source  locations 
These  standing  waves  propagate  across  H  field  lines  via  fast¬ 
mode  hydromagnetic  waves  Outside  the  surge  region  the 
waves  are  rapidly  damped  in  the  regions  of  lower  ionospheric 


conductivity.  This  model  therefore  provides  a  possible  energy 
path  that  Pi  2  pulsations  could  follow 

4  The  present  model  is  consistent  with  the  results  of  paper 
I  [Rothwell  et  al..  1984]  in  that  the  electron  precipitation 
fluxes  required  for  surge  propagation  are  also  sufficiently  high 
to  damp  the  excited  Pi  2  pulsations  The  flux  levels  predicted 
by  Fridman  and  Lemaire  p980]  also  ensure  the  presence  of 
damped  Pi  2  pulsations  in  the  model  given  here 

The  present  work  does  not  exclude  the  possibility  of  ad¬ 
ditional  magnetospheric  sources  of  Pi  2  pulsations.  We 
assume,  however,  that  these  Pi  2's  are  easily  reflected  by  the 
high-conductivity  region  and  that  the  ground-based  magne¬ 
tometers  most  efficiently  respond  to  the  ionospheric  source 
presented  here 
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